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This is a statistical analysis of the oceanographic time series measured across Fram Strait at a 
latitude of 78°50'N. Fram Strait is the deepest passage between the Arctic Ocean and the North 
Atlantic. There are up to 16 mooring lines with instruments at different depths measuring water 
temperature and velocity. These variables vary on different time scales and the challenge is to 
distinguish different spatial flow regimes. 

For Fram Strait, a temperature criterion is traditionally applied to identify water-masses, i.e. water 
volumes of similar origin. Interpolation leads to a vertical latitudinal 2D cross-section from which 
a scalar - the hypothetical area of waters within a certain temperature interval - can be extracted. 
The scalar is combined with a similar interpolation of the velocities to approximate the volume flows 
through the gateway. This approach is not only numerically expensive but also incorporates many 
assumptions. The present study suggest a new network-based approach to discriminate between 
flow regimes without the need to introduce artificial data through interpolation. The new approach 
not only reproduces the known flow patterns, but also reveals topographical features which are not 
captured by a standard water-masses analysis. 



Network theory offers a appealing framework to study 
local and non-local relationships of spatio-temporal vari- 
ables. Global climate networks are reconstructed from 
fields of observables such as ocean temperature at sea 
level 0. Network analysis of the temperature field un- 
covered wave-like structures of energy flow related to sur- 
face ocean currents. Having an unique set of highly 
resolved oceanographic time series in hand, network 
analysis is here used to gain insight into the flow pat- 
terns at Fram Strait. The data set is the first network 
which "goes into the deep" of the ocean, i.e. where the 
network is laterally extended. 

Fram Strait is approx. 300 km wide in the deep part, 
considering the (shallow) continental shelf leads to a 
width of 500 km. Therefore, the current study is one 
of the first regional climate networks [5], The data ba- 
sis here is an experimental measurement series instead 
of reanalysis data. The advantage is that every feature 
found by the network analysis can directly be judged by 
the experienced oceanographer. A yearly verification of 
the inter- and extrapolation is offered by a spatially high 
resolved hydrographic survey. The survey is carried dur- 
ing each cruise for deploying or recovering the moored 
instruments. It delivers a quasi-instantaneous snapshot 
of the water properties across Fram Strait (Fig. [2] shows 
the situation in 2011). In the following, the under- water 
buoys will be referred as moorings. 



I. THE ARCTIC OCEAN 

The Arctic Ocean is subject to changes e.g. as shrink- 
ing sea ice extent and the additional thinning of the ice. 
Warm impulses of water from the Atlantic Ocean are 
circulating in the Arctic basins. Water from the Pa- 
cific Ocean transports internal energy ("oceanic heat") 
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FIG. 1: Map of the mooring positions in Fram Strait. The 
meridional distance between 78.8°N and 79°N is around 20 km 
which is the same as the zonal spacing. For our analysis only 
the measurements from 2003 and afterward (o) are taken. 



through Bering Strait into the Arctic Ocean. Changes 
of sea ice extent, thickness and volume are observed by 
remote sensing and by drifting buoys with measuring 
equipment inside [3]. Monitoring the state of the water 
below the ice is technically difficult and requires inter- 
national efforts in terms of oceanographic surveys with 
ice-breakers or multiyear observatories deployed at the 
seafloor. 

The approach to gain synoptic (i.e. large-scale) insights 
about the Arctic Ocean exchanges is to monitor the in- 
and out-flow of water through Fram Strait. Although 
climatological scales are defined over periods of 30 years 
and more according to the World Meteorological Organi- 
zation (WMO), it seems possible to set up a preliminary 
energy- and volume balances: There are only four key 
gateways to the Arctic Ocean where water and sea-ice 
leaves or enters [4j . The Fram Strait is the only deep-sea 
opening of the Arctic Ocean to the other oceans, all other 
gateways are shallower and/or narrower, see Tab.! 
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FIG. 2: Temperature field of the upper 800 m of Fram Strait at a latitude of 78 °50' N. The data was acquired during a cruise 
with RV Polastern between June 24th and July 11th, 2011. White areas correspond to sea floor. The bold perpendicular lines 
symbolize the longitudes of the mooring lines, a degree longitude corresponds to roughly 20 km here. 



II. THE FRAM STRAIT TIME SERIES 

The two main currents of Fram Strait are a northward 
inflow of water from the Atlantic (next to Spitsbergen) on 
the one hand and a southward outflow of polar freshwater 
(next to Greenland) on the other hand. The East Green- 
land Current (EGC) can be recognized as (dark blue) 
cold patch above 200 m at the left side of Fig.[2j The wa- 
ter there moves southwards (out of the paper plane) and 
carries sea ice out of the Arctic Ocean. The West Spits- 
bergen Current (WSC) is the patch with temperatures 
above 2° C (red, yellow, light-green) on the right side of 
Fig.[2] east of 5° E. The WSC is a branch of the North 
Atlantic Current which is the northern extension of the 
Gulf Stream. The WSC is the warmest water mass en- 
tering the Arctic Ocean. Since our measurements started 
in 1997, the temperature of the Atlantic water increased 
with a rate of around 1° C within 10 years [3]. 

Motivated by the importance of the Fram Strait for the 
Arctic climate, the Alfred- Wegener-Institute for Polar- 
and Marine Research (AWI) maintains a transect of 
moorings in collaboration with the Norwegian Polar In- 
stitute at 78°50'N. The transect ends at the shelf breaks 
at 6°52W and 8°40E respectively, see Fig. [2] and Fig. [TJ 

Instruments (from the companies Aanderaa and Fal- 



Arctic gateway 


width 


max. depth 


net volume flux 


Barents Sea Opening 


350 km 


500 m 


2.0±2.9Sv 


Fram Strait 


300 km 


2 600 m 


-2.0±2.7Sv 


Bering Strait 


86 km 


50 m 


0.8±0.2Sv 


Nares Strait 


40 km 


220 m 


0.57 ± 0.3 Sv 


Lancaster Sound 


40 km 


125 m 


-0.7Sv 


Cardigan Strait 


8 km 


180 m 


-0.3 Sv 



TABLE I: Gateways to the Arctic Ocean, the last three are 
the largest of the Canadian Archipelago. 1 Sverdrup (Sv) = 
10 6 m 3 /s. Emptying Lake Constance (German: "Bodensee") 
within a day corresponds to 0.56 Sv. Estimates from [4], pos- 
itive values: inflow into the Arctic Ocean. 



mouth Scientific Instruments) provide point measure- 
ments of horizontal velocities, flow direction, tempera- 
ture and salinity. However, the salinity measurements 
of the current meters are not trustworthy since the con- 
ductivity cell is not pumped. This causes much slower 
response times than what is needed to cope with the 
high flow speeds in Fram Strait. Besides, the sensors 
fail frequently due to growing bio-films. Therefore, many 
gaps within the salinity time series exist. The addition- 
ally available data from a few upward-looking Acoustic 
Doppler current profilers (ADCP) near the surface are 
also neglected as their time-series are far shorter. 

Data is recorded at 2 to 5 different depth levels every 
hour. The lowest level usually is located approximately 
10 m above the sea-floor. If applicable, there are instru- 
ments at 1500 m depth, around 800 m and 300 m depth. 
The uppermost sensor is located around 60 m below the 
sea surface to avoid losses of instruments by drifting ice 
keels. Recently, fish trawling became more frequent with 
the retreating ice edge. The fishing nets are a threat for 
our instruments. 

In-situ data by observational oceanography are assim- 
ilated by ocean modelers. A major product for them are 
time series of de-tided daily averages. The effect of tides 
are excluded by filtering out semi-diurnal and diurnal 
tidal constituents. The analysis is based on the daily 
means, accepting that any additional post-processing 
of the measurements (e.g. depth-correction to due flow- 
induced dives or manual spike-removal) happens inside a 
black box. For our purposes these are only minor changes 
which do not change the main statements. 



III. CONSTRUCTING A CORRELATION 
NETWORK 

The method suggested here should work for any set of 
spatio-temporal time series. In Fram Strait, the network 
of point measurements has 66 nodes with three times 
series per node x: the temperature T x , the meridional 
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FIG. 3: Temperature correlations coefficients between nearest and 2nd nearest neighbors of all timeseries since 2003 for daily 
means (top) and with a 7-day filter (bottom). The bullets represent the average instrument position (depth and longitude), for 
this reason the mooring lines are not straight as sketched in Figj2] 
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velocity v x and the zonal velocity u x at a given posi- 
tion x G {1, . . . , 66}. Each x can be mapped on a lon- 
gitudinal position and a depth level. Every time series 
is constructed by merging snippets of the maintenance 
period. 

The Spearman correlation coefficient p between any 
nearest or second-nearest neighbors is calculated for the 
complete time series since 2003. Working with daily av- 
erages, only correlations between observables at the same 
time instance are considered. Pairwise complete observ- 
ables are taken into account, i.e. whenever both time se- 
ries have a valid value at a given point in time. Looking at 
all different pairs of correlation is a standard consistency 
check and therefore 66 2 — 4365 correlation coefficients 
have to be evaluated. New is the concise and intuitive 
representation as a network of (2nd) nearest neighbors 
which add up to only 205 correlation coefficients. 

One expects high correlations even between non- 
neighboring pairs of point measurements because the dy- 
namics of the water is driven by the atmosphere and 
ocean upstream. The new method is to use a multi-scale 
approach which does not look on the value at the cor- 
relation itself, but at how the correlation increases with 
increasing filtering size. Before, it has to be noted that 
correlations between neighboring instruments are high 
in the vertical and significantly lower in the horizontal 
direction, simply because of the different resolutions in 
both directions. Within a 2 km mooring line there are 
up to five instruments while there is only a mooring line 
roughly every 20 km. The higher vertical correlations 
(between velocities within same mooring line) are usu- 
ally interpreted as sign of a prevailing barotropic struc- 
ture of the flow, i.e. the flow is mainly directed by the 
topography of the ground and does not depend on the 
depth. 

Moving towards longer filtering times, the zonal corre- 
lation within the (sub) surface widens towards the West. 
This is expected when the upper ocean is well-mixed by 
the atmosphere. In particular, where the ocean surface 
is only rarely covered by sea ice, mixing takes place. The 
strong correlations at the bottom between 2 and 4° East 
are most probably induced by the topographic steering of 
the Knippovich Ridge whose Northern foothill mountains 
are located at these longitudes, see Fig. [I] Besides, there 
are vertical correlations at the correct depth and longi- 
tudinal position where the Atlantic Water is expected 
to recirculate. The recirculation manifests itself in form 
of coherent southward flow. A hint is the warm patch 
(T > 4°C) between ±3°E in Fig.[5] (REC). The long- 
term averaged velocity field (not shown here) supports 
this assumption. 

The effect of the filtering time on the pattern is small 
for the correlation network constructed from the water 
velocities. Therefore, without loss of generality a 10-days 
boxcar filter is applied in Fig. [4] 



IV. MULTISCALE CORRELATION ANALYSIS. 

In oceanography, a water mass is "an identifiable water 
volume with a common history which may be traced back 
to some source area" . Traditionally, water masses are de- 
fined by (potential) temperature and salinity, although 
additional properties may also taken into account. Ide- 
ally, these are conservative, in the sense that they are 
only modified by mixing, and neither biological activity, 
nor chemical degradation influences the parameter. In 
Fram Strait, no reliable salinity measurements for the 
whole time series exist. Thus, a pure temperature cri- 
terion has been introduced to identify different water 
masses. Whenever the temperature is above a threshold 
of 1°C the water is assumed to originate in the Atlantic 
[BJ. Here, a statistical method to identify water masses 
is presented. This approach can be also helpful for other 
networks of time series. 

The correlation pattern between the 2nd and 3rd moor- 
ing line from the East is independent of the variable, see 
Fig. [3] and [4] It is exactly the location of the central core 
of the West Spitsbergen Current (WSC) which is con- 
sistent with the temperature-criterion for identifying the 
water masses, compare Fig.[2j 

It is even more interesting to identify regions where 
correlations are lower and the absence of links persists 
when changing the filtering time. All correlations grow 
slowly with increasing filtering intervals since the filter- 
ing operation introduces correlations. However, the effect 
of growing filter window also provides oceanographic in- 
formation. For instance the missing vertical links in in 
the Western part of the Strait (for all variables, but more 
pronounced for the temperature-correlation network) are 
a clear sign for the decoupling of the ocean from the at- 
mosphere which occurs due to frequent sea ice cover in 
this region. For the water velocities in Fig.|4j the oth- 
erwise prominent vertical correlations disappear between 
the 2nd and 3rd layer of instruments. Exactly in those 
regions the two boundary currents tend to have their ver- 
tical limits. 



V. CONCLUSION 

Traditional methods to classify flow regimes rely on 
numerical expensive interpolation schemes. The draw- 
back is that values (for temperature and water velocity) 
have to produced for the complete vertical field across the 
gateway. The interpolation is based on sparse measure- 
ments every 20 km which necessarily comes along with 
large error bars for transport estimates, see pQ for a reli- 
able method based on finite elements. 

A new climate network has been introduced which al- 
lows to identify water masses and other spatially con- 
sistent flow regimes using simple maths. The nodes are 
the instrument positions, given as longitude-depth coor- 
dinates. Links are established if the correlation between 
the time series between two nodes is bigger than a certain 
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FIG. 4: Correlation graph for meridional (top) and zonal (bottom) velocities with a 10-day filter for the complete time series 
since 2003. 



threshold. 

The new approach allows to quantify regions of in- 
terest without any prior knowledge about oceanographic 
phenomena. Where links are absent within the corre- 
lation network, increasing the spatial resolution of the 
transect may help to gain precision. One crucial region is 
the longitudinal resolution at 6±1° E. The new method 
quantifies the common knowledge that the error of the 
flux estimates actually originates in the uncertainty of 
the western boundary of the West Spitsbergen Current 
(WSC). 



VI. DATA SOURCES 

Eventually the complete data set will be freely avail- 
able, the following table summarizes the state in June 
2012. 



snippet begin end °W °N doi:10.1594/ 
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78 


832 


pangaea 


778858 


F10-7 


2004- 


08- 


24 


2005 


08- 


29 


-2.001 


78 


831 


pangaea 


778859 


F10-9 


2006- 


09- 


09 


2008 


07- 


21 


-2.050 


78 


821 


pangaea 


777567 


F10-10 


2008- 


07- 


21 


2010 


07- 


20 


-2.115 


78 


828 


pangaea 


777566 


F15-2 


2003- 


09- 


28 


2004 


08- 


23 


1.611 


78 


833 


pangaea 


778860 


F15-3 


2004- 


08- 


23 


2005 


08- 


29 


1.610 


78 


833 


pangaea 


778861 
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08- 


30 


2006 


08- 


30 


1.610 


78 


833 
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08- 


30 
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09- 


24 


1.609 


78 


833 


pangaea 


777568 
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2007- 


09- 


24 


2008 


07- 


19 
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78 


833 


pangaea 


777569 
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2008- 


07- 


18 


2010 


07- 


17 


1.599 


78 


834 


pangaea 


777570 


F16-2 


2003- 


09- 


29 


2004 


08- 


22 


0.401 


78 


835 


pangaea 


778863 


F16-3 


2004- 


08- 


22 


2005 


08- 


29 


0.397 


78 


834 


pangaea 


778864 


F16-4 


2005- 


08- 


30 


2006 


08- 


30 


0.401 


78 


835 


pangaea 


778865 


F16-5 


2006- 


08- 


31 


2007 


09- 


13 


0.401 


78 


835 


pangaea 


777571 


F16-6 
F16-7 


2007- 
2008- 


09- 
07- 


13 

20 


2008 
2010 


07- 
07- 


19 
17 


0.540 
0.390 


78 
78 


832 
827 


pangaea 


777572 


pangaea 


777573 



The first column is the identifier of the snippet, the 
first number encodes the position, the number after the 
dash counts the number of deployments. 
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